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In  this  work  OH  observations  of  the  middle  atmosphere  per¬ 
formed  by  the  improved  2.5 THz  heterodyne  spectrometer  THOMAS 
(Tera  Hertz  OH  Measurement  Airborne  Sounder)  during  the  1997 
THOMAS/MAHRSI  (Middle  Atmosphere  High  Resolution  Spectrograph 
Investigation)  campaign  are  presented.  Residts  of  the  THOMAS  OH 
measurements  arc  compared  to  simultaneous  MAHRSI  OH  observations  and 
to  photochemical  model  calculations  using  both  standard  HO*  chemistry 
and  a  recently  proposed  change  in  IIOx  chemistry.  The  10  km  90  km 
OH  column  densities  measured  by  THOMAS  and  MAHRSI  show  excellent 
agreement,  the  50  km  -  90km  OH  column  densities  agree  well  within  the 
measurement  uncertainties  throughout  the  diurnal  cycle.  The  comparison 
of  the  observed  and  modeled  OH  column  densities  in  the  upper  stratosphere 
and  mesosphere  suggests  that  neither  model  is  capable  of  simultaneously 
reproducing  the  measurements  in  both  regions. 


1.  INTRODUCTION 

Chemical  reactions  involving  the  HO*  family  (H, 
OH,  and  H02)  play  a  significant  role  in  middle  atinos- 
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pheric  ozone  chemistry.  However,  recent  measurements 
of  HO*  in  the  upper  stratosphere  and  mesosphere  present 
an  inconsistent  picture,  indicating  either  measurement 
errors  or  deficiencies  in  our  understanding  of  atmos¬ 
pheric  odd  hydrogen.  Up  to  now,  measurements  of 
middle  atmospheric  OH  could  not  be  compared  di¬ 
rectly  because  they  were  carried  out  at  different  loca¬ 
tions,  at  different  times,  or  covered  altitude  regions 
that  do  not  overlap.  For  instance,  retrievals  from  the 
1994  space  borne  observations  of  mesospheric  OH  by 
MAHRSI  [ Conway  et  al.,  1996,  Summers  et  al.,  1997, 
Conway  et  al.,  1999)  extend  down  to  50km,  whereas 
t  he  measurement  of  stratospheric  OH  profiles  by  a  bal¬ 
loon  borne  instrument  (FIRS  2)  covers  altitudes  below 
50km  [Jacks  et  al..  1998].  For  the  balloon  borne  meas- 
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urement,  the  reported  OH  concentrations  at  50  km  are 
in  agreement  with  the  standard  chemistry  prediction, 
while  the  1994  MAHRSI  results  were  found  to  In*  .30  91 
-  40  %  lower  than  the  calculated  model  values  at  50  km 
\WM(),  1998],  To  resolve  these  problems,  coordinated 
atmospheric  measurements  with  different  measurement 
techniques  and  overlapping  altitude  regions  are  desir¬ 
able  in  order  to  rule  out  potential  unknown  systematic 
errors. 

OH  measurements  covering  the  middle  and  upper 
stratosphere  and  the  mesosphere  wpre  performed  in  Au¬ 
gust  1997  with  the  improved  airborne  2.5  TH*  hetero¬ 
dyne  spectrometer  THOMAS.  These  measurements  were 
coordinated  in  space  and  time  with  measurements  of 
mesospheric  and  upper  stratospheric  OH  from  the 
second  MAHRSI  flight.  In  the  following,  we  give  a 
brief  overview  of  the  THOMAS  instrument,  the  in¬ 
version  technique  and  the  method  applied  to  compare 
THOMAS  observations  to  MAHRSI  results.  Finally, 
the  results  of  the  THOMAS  OH  measurements  are  pre¬ 
sented  and  compared  to  MAHRSI  measurements  as  well 
as  to  photochemical  model  calculations. 

2.  EXPERIMENT 

2.1.  Instrument 

The  THOMAS  instrument  is  a  further  development 
of  a  2.5  TH i  heterodyne  spectrometer  for  extraterres¬ 
trial  observations  that  was  built  at  the  Max  Planck 
Institute  for  Radio  Astronomy  in  Bonn,  Germany  [/loser. 
1991).  At  DLR  (German  Aerospace  Center.  Ober- 
pfaffenhofen),  the  instrument  design  was  modified  to 
allow  the  observation  of  OH  thermal  emission  lines 
in  earth's  atmosphere  from  the  DLR  research  aircraft 
FALCON,  using  an  up-looking  geometry.  The  oper¬ 
ation  of  THOMAS  above  the  tropopause  is  desirable 
because  of  the  significant  signal  absorption  due  to  wa¬ 
ter  vapor  in  the  far  infrared  at  lower  altitudes.  The 
first  measurement  flights  were  performer!  in  1994/95, 
proving  the  airworthiness  of  the  instrument  [  Titz  ct  of., 
1995a;  Titz  et  at ,  1995b].  Between  these  flights  and  the 
1997  THOMAS  measurements,  the  spectrometer  was 
significantly  improved  and  a  new  aircraft  window  was 
built  [fc'np/ert  et  al..  1999a:  Englert  et  of.,  1999b;  Eng- 
lert,  1999). 

The  instrumental  modifications  increased  the  signal 
to  noise  ratio  by  more  than  a  factor  of  five,  while  im¬ 
proving  the  spectral  resolution  by  about  a  factor  of 
two.  Figure  1  show’s  a  typical  spectrum  measured  by 
THOMAS  after  the  improvements.  A  strong  water 
vapor  line  and  the  OH  triplet  (two  strong  and  one  weak 
line  near  83.869cm-1)  can  clearly  be  recognized. 
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Figure  1.  Calibrated  spectrum  measured  during  the 
THOMAS/ MAHRSI  campaign  in  August  1997 

The  photometric  accuracy  of  THOMAS  was  investi¬ 
gated  using  laboratory  measurement  s  of  methanol  emis¬ 
sion  lines  and  was  included  in  the  error  budget. 

2.2  Measurement  Geometry  ami  Time  Coverage. 

The  measurement  locations  and  timing  were 
selected  to  achieve  nearly  zero  miss  distance  and  zero 
miss  time  with  respect  to  the  MAHRSI  observations  for 
the  upper  stratosphere  and  mesosphere.  Due  to  the  dif¬ 
ferent  viewing  geometries  of  the  instruments  (THOMAS: 
up-looking,  70°  85°  from  zenith,  MAHRSI:  limb  view¬ 

ing,  tangent  heights  between  about  35  km  and  9t)km)  a 
perfect  match  of  the  probed  air  volume  is  not  possible. 
However,  since  t  he  upper  stratospheric  and  mesospherir 
OH  column  amount  is  dominated  by  Oil  between  41)  km 
and  60  km.  the  FALCON  flight  route  was  planner]  in  a 
way  that  the  air  volume  observed  by  THOMAS  at  50  km 
coincided  with  the  tangent  point  track  of  MAHRSI. 
whenever  possible.  Plate  1  shows  the  measurement 
geometry  for  August  15,  1997  The  two  lines  indi¬ 
cating  the  layer  entrance  and  exit  points  of  a  virtual 
MAHRSI  30  km  tangent  height  line  of  sight  (layer  thick¬ 
ness:  5  km)  were  added  to  the  plate  to  illustrate  the  ex¬ 
tent  of  the  air  volume  measured  by  MAHRSI  Minimal 
miss  distance  and  time  was  achieved  for  ten  MAHRSI 
overflights  on  4  days  between  August  10,  and  August 
15,  1997.  The  coordinated  measurements  range  over  a 
local  solar  time  (LST)  interval  between  about  7 am  and 
5pm.  An  additional  sunset  flight  was  performed  with 
THOMAS  to  ejetend  the  time  coverage  to  about  8  pm 

LST. 

3.  DATA  ANALYSIS 

31.  Inversion 

A  best  fit  of  a  radiative  transfer  calculation  to  the 
measured  spectrum  is  used  to  determine  OH  concen¬ 
tration  profiles  from  the  up  looking  thermal  emission 
spectra.  The  information  about  the  vertical  distribu¬ 
tion  of  OH  is  contained  in  the  pressure  dependent  line 
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width  (OH  line  parameters  were  adopted  from  Park  et 
al.  (1906]  (pressure  broadening  and  temperature  expo¬ 
nent)  and  the  HITRAN  data  base  (line  position  and 
strength)).  The  inversion  of  concentration  profiles  from 
atmospheric  emission  data  is  generally  known  to  be  an 
ill- poser  1  problem.  Regularization ,  i.e.  the  introduction 
of  additional  constraints,  is  necessary  to  get  a  numen 
cally  stable  and  physically  meaningful  solution. 

The  THOMAS  spectra  have  been  processed  using 
Phillips-Tikhonov  regularization  (PTR)  as  developed 
by  Hansen  [1992]  and  first  applied  to  atmospheric  re¬ 
trieval  by  Schtmpf  and  Schrtier  [1997].  In  contrast 
to  optimal  estimation  [Rodger*,  1976]  PTR  does  not 
require  specification  of  an  a  priori  profile  and  its  co- 
variance  matrix-  Instead,  PTR  stabilizes  the  solution 
of  the  least  squares  problem  by  adding  a  smoothness 
constraint,  and  Hansen's  L-curve  criterion  is  used  to 
automatically  determine  the  optimal  smoothing,  allow¬ 
ing  the  extraction  of  as  much  information  as  possible 
from  the  measurement  | Sehtmpf.  1999,  Englert  et  al. 
manuscript  in  preparation ], 

In  addition  to  the  OH  concentration  profile,  the  in¬ 
version  algorithm  also  calculates  the  averaging  kernel 
matrix.  In  theory,  the  inverted  spectrum  is  the  product 
of  the  averaging  kernel  matrix  and  the  true  concentra¬ 
tion  profile,  which  means  that  the  averaging  kernel  ma¬ 
trix  rows  determine  the  sensitivity  of  the  ret  rieved  OH 
concentration  at  a  particular  altitude  to  the  OH  con¬ 
centration  at  other  altitudes.  Due  to  the  smoothing  of 
the  inverted  profile  by  the  averaging  kernel,  it  is  com 
monly  used  to  characterize  the  altitude  resolution  of  the 
retrieved  profile.  The  middle  panel  of  Plate  2  shows  typ¬ 
ical  THOMAS  averaging  kernel  rows,  illustrating  how 
the  resulting  profile  is  smoothed  on  the  altitude  grid  of 
the  inversion.  In  the  middle  and  upper  stratosphere, 
the  averaging  kernel  rows  are  strongly  peaked  and  the 
typical  full  width  at  half  maximum  of  the  rows  is  around 
25  km.  In  the  mesosphere  however,  the  averaging  ker¬ 
nel  rows  are  not  strongly  peaked,  showing  that  in  this 
region  no  hpight  information  is  available  from  the  spec¬ 
trum.  This  behavior  is  a  consequence  of  the  fact  that  in 
the  mesosphere  the  line  shape  of  OH  thermal  emission 
is  no  longer  dominated  by  pressure  broadening,  but  by 
temperature  broadening,  which  results  in  a  nearly  alti¬ 
tude  independent  line  shape  in  the  mesosphere. 

3.2.  Determination  of  Weighted  Column  Densities 

Due  to  the  substantial  smoothing  of  the  Oil  profiles 
inverted  from  THOMAS  data,  a  direct  comparison  to 
high  vertical  resolution  MAHRS1  OH  profiles  or  photo¬ 
chemical  model  results  is  not  possible.  Multiplication 
of  the  high  vertical  resolution  profiles  with  a  THOMAS 


averaging  kernel  can  be  used  to  match  the  vertical  re¬ 
solution  of  the  profiles  to  be  compared.  The  compar¬ 
isons  presented  here  utilize  weighted  column  densities , 
defined  as  the  atmospheric  Oil  profile  smoothed  by  a 
THOMAS  averaging  kernel  and  integrated  over  altitude 
as  illustrated  in  Plate  2.  The  left  panel  of  Plate  2  shows 
two  photochemical  model  OH  profiles  representing  the 
standard  odd  hydrogen  chemistry  \Dt\1ore  et  at.,  1997] 
(henceforth  referred  to  as  Model  A)  and  revised  odd  hy¬ 
drogen  chemistry  (henceforth  referred  to  as  Model  B) 
after  Summers  et  al.  [1997],  III  the  right  panel  the 
same  profiles  are  shown  after  smoothing  with  the  aver¬ 
aging  kernel  winch  is  depicted  in  the  middle  panel  In 
addition  to  the  model  profiles,  a  THOMAS  measured 
profile  is  included  in  the  right  pm  id  of  Plate  2.  The 
shaded  area  illustrates  the  40  km  90  km  weighted  OH 
column  density.  The  smoothing  of  the  OH  profiles  prior 
to  the  integration  makes  the  wetghted  column  densities 
dependent  on  OH  below  the  lower  limit  of  integration. 
For  the  50  km  to  90  km  weighted  OH  column,  the  sen¬ 
sitivity  to  OH  at  lower  altitudes  is  about.  20%  whereas 
it.  is  about  15%  for  the  40  km  to  90  km  wciglitisl  OH 
column. 

The  different  altitude  coverage  of  THOMAS  and 
MAHRSI  further  complicates  the  comparison  of  the  re¬ 
sults  of  the  two  experiments.  OH  profiles  inverted  from 
individual  MAHRSI  limb  scans,  that  are  used  for  the 
minimum  miss  time  and  distance  comparisons,  cover 
the  altitude  region  between  about  50  km  and  90  km  In 
order  to  smooth  the  single  scan  MAHRSI  profiles 
using  a  THOMAS  averaging  kernel,  the  MAHRSI  pro¬ 
files  have  to  be  extended  into  the  stratosphere.  The 
extension  into  tin*  upper  stratosphere  was  performed 
based  on  averaged  MAHRSI  limb  scan  data  [ Conway 
et  al.,  2000]  that  allowed  a  reliable  inversion  of  OH 
concentrations  down  to  about  38  km  Below  38  km 
the  extension  was  made  using  standard  photochem¬ 
istry  results  (Model  A)  that  show  reasonable  agreement 
with  MAllIlSi  results  at  38km  [Conway  et  at.,  2000 
Using  this  technique,  50km  to  90  km  weighted  OH  col¬ 
umn  densities  have  been  derived  from  single  MAHRSI 
limb  scans.  40  km  to  90  km  weighted  column  densities 
have  heen  derived  from  averaged  MAHRSI  limb  scans 
extended  by  Model  A  results  below  38  km. 

3.3.  Error  Assessment 

The  assessment  of  possible  errors  in  the  THOMAS 
weighted  OH  column  amounts  was  performed  using  a 
typical  noise  contaminated  synthetic  spectrum.  The 
synthetic  spectrum  has  been  inverted  including  pertur 
ballons  due  to  assumed  error  sources,  in  order  to  esti¬ 
mate  the  effects  on  the  weighted  OH  column  density. 
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Plate  1.  Measurement  geometry  of  August  15,  1997.  The  ’predicted’  data  is  the  result  from  initial 
calculations. 


Plate  2.  Left  panel:  Model  OH  profiles  using  the  standard  (Model  A)  and  a  revised  chemistry  (Model  B). 
Center  panel:  typical  THOMAS  averaging  kernel  representing  the  mapping  of  the  high  resolution  model 
profiles  to  the  low  resolution  domain  of  THOMAS  retrieval  results.  The  rows  correspond  to  fifteen 
altitude  layers  of  equal  spacing  between  15  km  and  85  km.  Right  panel:  Both  model  profiles,  smoothed 
with  the  averaging  kernel  so  that  they  can  be  compared  directly  to  the  corresponding  profile  retrieved 
from  a  THOMAS  measurement  (LST  «  13:00).  (Error  bars  illustrate  the  statistical  error  only.)  The 
shaded  area  represents  the  weighted  column  density  between  40  km  and  90  km. 
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Figure  2.  Comparison  of  40  km  —  90  km  weighted  OH 
column  densities  derived  from  THOMAS  observations, 
MAHRSI  results  from  averaged  limb  scans,  and  two  pho¬ 
tochemical  model  calculations  The  MAHRSI  profiles  from 
averaged  limb  scans  cover  the  38  km  -  90  km  altitude  region. 
Mudel  A  uses  standard  odd  hydrogen  chemistry  [PeMnrr.  e.t 
ai,  1997],  model  0  uses  a  50%  reduction  in  the  rate  coef¬ 
ficient  fui  O-eHOj  — >  Oll-eOj  a*  discussed  in  Surameo  e< 
of.,  [1997]. 

For  the  1997  THOMAS  measurements,  eighteen  sys¬ 
tematic  error  sources  are  investigated  this  way  and  the 
total  systematic  error  for  the  40  km  -  90  km  as  well 
as  the  50  km  90  km  weighted  Oil  column  density  is 
found  to  be  on  the  order  of  ±15%.  The  error  sources 
regarded  here  include  instrumental  errors,  uncertainties 
in  atmospheric  pressure  and  temperature,  variations  in 
interfering  gas  concentrations  and  uncertainties  in  spec¬ 
tral  line  and  continuum  data.  The  statistical  error  of 
the  inversion  is  typically  in  the  range  of  ±7%. 

4.  RESULTS 

Figure  2  shows  40  km  90  km  weighted  column  den¬ 
sities  as  a  function  of  I. ST  that  were  derived  from 
THOMAS  data,  averaged  MAHRSI  limb  scans  and  from 
photochemical  model  results.  Model  A  uses  standard 
odd  hydrogen  chemistry  whereas  Model  B  incorporates 
a  50%  decrease  in  the  key  OH/HOj  partitioning  reac¬ 
tion  (O+HOj  OH  •+■  O2)  as  discussed  by  Summers 
et  al..  [1997],  The  comparison  of  the  40  km  to  90  km 
weighted  OH  column  densities  shows  excellent  agree¬ 
ment  of  THOMAS  and  MAHRSI  observations.  The 
upper  stratospheric  and  mesospheric  OH  column  abun¬ 


dances  observed  by  both  instruments  are  about  15% 
lower  than  the  model  prediction  using  the  standard  odd 
hydrogen  chemistry  [Model  A).  The  model  results  using 
the  revised  chemistry  are  about.  10  %  lower  than  the  ex¬ 
perimental  values. 

Figure  3  shows  50  km  90  km  weighted  OH  column 
densities  derived  from  THOMAS  observations, 
MAHRSI  single  scan  observations  coinciding  in  space 
and  time  with  the  THOMAS  observations  and  photo¬ 
chemical  model  resu-ts  as  a  function  of  LST.  The  50 km 
-  90  km  weighted  OH  column  densities  measured  by 
THOMAS  are  abou:  15%  lower  than  Model  A  results 
and  about  15%  greater  than  predicted  by  Model  B. 
The  values  derived  from  coincident  MAHRSI  single  limb 
scans  show  agreement  with  THOMAS  to  better  than 
about  10%,  but  are  generally  lower  than  the  THOMAS 
observations. 

The  overall  shape  of  the  diurnal  variation  of  the  OH 
weighted  column  densities  depicted  in  Figures  2  and  3  is 
in  agreement  with  mode]  predictions.  This  shows,  that 
neither  local  variations  in  OH,  nor  changes  between  the 
days  of  observation  have  been  significant  regarding  the 
statistical  error  of  the  presented  experimental  data. 

The  error  bars  depicted  in  Figures  2  and  3  represent 
the  statistical  errors  due  to  measurement  noise  only.  As 
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Figure  3.  Comparison  of  50  km  -  90  km  weighted  OH  col¬ 
umn  densities,  derived  from  THOMAS  observations,  single 
limb  scan  MAHRSI  results,  and  two  photochemical  model 
calculations.  The  single  limb  scan  MAHRSI  profiles  cov¬ 
er  the  50km  90km  altitude  region  and  coincide  in  space 
and  time  with  the  respective  THOMAS  measurements.  The 
denotation  of  the  model  results  is  analogous  to  Figure  2. 
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discussed  in  Section  3.3,  the  maximum  systematic  er¬ 
ror  for  the  OH  column  amounts  observed  by  THOMAS 
is  approximately  15%.  The  total  systematic  errors  of 
the  MAHRSI  measurements  are  estimated  to  be  1 12.0/ 
-15  5%  for  the  50  km  90  km  column  and  +  !9.0/-22.5% 
for  the  -10  km  90  km  column  density  (see  also  (Con tray 
et  al.,  1999]) 

5.  CONCLUSIONS 

The  results  of  this  first  joint  upper  stratospheric  and 
mesospheric  OH  measurement  campaign,  using 
two  entirely  different  measurement  techniques  (different 
geometry  of  observation  and  different  spectral  regions), 
show  that  measurements  of  mesospheric  Oil  are  gen¬ 
erally  lower  than  expected  from  standard  photodicin- 
ical  theory.  This  conclusion  brings  into  question  our 
understanding  of  the  upper  stratospheric  ll()r  chem¬ 
istry,  since  the  same  HO,  chemical  reactions  dominate 
OH/HO2  partitioning  in  both  regions.  In  view  of  the 
latest  balloon  observations  by  Jacks  et  al  (1998),  long 
term  column  measurements  by  Burnett  ct  til.  [1998]  and 
the  1994/97  MAIIRSI  results,  it  seems  unlikely  that  a 
rate  change  of  a  single  reaction  rate  in  H0r  chemistry 
will  lead  to  agreement  between  HO,  models  ai.d  meas¬ 
urements  throughout  the  middle  atmosphere. 
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